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Motivated by recent experiments [T] [2], we study quasi 2D ferromagnetic condensates with various 
aspect ratios. We find that in zero magnetic field, dipolar energy generates a local energy minimum 
with all the spins lie in the 2D plane forming a row of circular spin textures with alternating 
orientation, corresponding to a packing of vortices of identical vorticity in different spin components. 
In a large magnetic field, the system can fall into a long lived dynamical state consisting of an array 
of elliptic and hyperbolic Mermin-Ho spin textures, while the true equilibrium is an uniaxial spin 
density wave with a single wave- vector along the magnetic field, and a wavelength similar to the 
characteristic length of the long lived vortex array state. 



The condensates of bosons with non-zero spins, known 
as spinor condensate, are remarkable superfluids. In ad- 
dition to broken gauge symmetry, they also have broken 
symmetries in spin space. The spin degrees of freedom 
lead to a variety ground states, which proliferates rapidly 
as the value of spin increases. Since different spin com- 
ponents can be mixed through spin rotation, there is 
considerable interplay between spin and gauge degrees 
of freedom, leading to a whole host of new macroscopic 
quantum phenomena. 

The simplest spinor condensates are those for spin- 
1 bosons, such as the F = 1 hyperfine states of ^^Na 
and ^^Rb. The ground state of ^^Na is a non-mangetic 
"polar" condensate whereas ^^Rb is a ferromagnetic 
condensate [3]. In the case of ferromagnetic condensates, 
they possess an additional "spin-gauge" symmetry which 
makes non-uniform spin textures behave like vortices [4^. 
The system can respond to external rotation through spin 
deformation. Any attempt to bend the spin will also gen- 
erate vorticity. 

The magnetic nature of spinor condensates naturally 
leads to the consideration of dipolar energy, which is in- 
trinsic to alkali atoms. Since dipole energy can generate 
non- uniform spin textures, it will generate vorticity. In- 
deed, Yi and Pu have shown that a ^^Rb condensate in 
a sufficiently flat cylindrical potential will form a circu- 
lar spin texture, which is a vortex of ferromagnetic con- 
densate with a polar core [5 . Recently, experiments at 
Berkeley have shown that a ^^Rb condensate with a he- 
lical texture can decay into a random spin textures [6]. 
By estimating the energy of the final state, the authors 
suggest that the phenomenon is caused by dipolar en- 
ergy. More recently, the Berkeley group has found that 
a pancake like condensate of ^^Rb can develop a texture 
with periodically modulated spin-spin correlation rotat- 
ing rapidly about an in-plane magnetic field [U [2]. They 
suggest that this effect is also due to dipolar energy. 

Dipolar effects are highly geometry dependent. In 
this paper, we would like to point out that some key 
features of quasi 2D ^^Rb condensate due to dipolar 
interactions. Much of what we discuss also apply to 



other ferromgnetic condensates. We shall consider an 
anisotropic trap with frequencies uOz » oOy > uOx^ as in 
ref. m HI E] • The condensate is then a thin anisotropic 
slab in the x?/-plane with Thomas- Fermi radii R^^Ry 
such that Rx/Ry = [ojy/u^j,) = A, (A - 10 In ref. [Q El [6]). 
In our discussions, we choose the normal to the conden- 
sate slab, z, to be the spin quantization axis for the con- 
densate wavefunction = ?/;o, V^-i), where the su- 
perscript "T" stands for transpose. We shall show that 

(1) In zero magnetic field, dipolar energy leads to a lo- 
cal energy minimum consisting of a row of circular spin 
textures with alternating spin orientations in the long 
direction x, with all the spins in the xy-plane. These 
textures are of the Yi-Pu type [5 , with a polar core and a 
size determined by the Thomas- Fermi radius in the short 
direction, Ry. This state amounts to an array of vortices 
in identical vorticity in the ipi component, with be- 
ing its time reversed partner. We have also found an 
analytic expression that well approximates this state. 

(2) In a large magnetic field B in the xy-plane, the spins 
rotate about B and experience a time averaged dipolar 
energy. We find that this energy is very flat in spin space 
around a class of textures which is an array of elliptical 
and hyperbolic Mermin-Ho vortices. These states are not 
local minima. They will eventually evolve to the true 
minimum, which is an a uniaxial spin density wave (or 
"stripe phase" for short) with a single wave-vector along 
B. Due to the flatness of the energy surface, the vortex 
lattice textures will be very long lived during dynamical 
evolutions. 

(A) Basic Structures: We first consider some spin 
textures relevant for later discussions. We shall write the 
condensate wavefunction as (x) = ^/n{x.)(^ (x) , where 
n(r) = l^/iP is the density, and C*" • ^ = 1. The spin 
field is given by S = ^* F^^^^^y = n(x)m(x) where F^j^ is 
the spin operator, and m = C/^F^^Czy- The general form 
of a ferromagnetic condensate is 

C = e '^ [u ,v2uv,v or — - — = —= (1) 
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FIG. 1: (la) and (lb) are elliptic planar textures with iden- 
tical vorticity in ipi and ip-i. (Ic) and (Id) are hyperbolic 
planar spin textures with identical vorticity. All these spin 
textures (la)-(ld) have a polar core, (le) and (If) are ellip- 
tic and hyperbolic Mermin-Ho textures. The cores are fully 
ferromagnetic. 



where u = e~*"/^cos/3/2, v = e*"/^sin/3/2. The spin is 
m = cos/3z + sin/?(cosax + sinay), and = 1. If the 
spin lies in the x?/-plane, then \u\ = \v\ = 1/a/2, (or 
IV^iP : IV^oP : IV^-iP = 1 • 2 : 1), C then reduces to 



^{e a/2, e*")/2. The following cases are of 



interest to us. (Below, (r, (/)) are polar coordinates.) 
(i) Elliptic planar texture: This corresponds to 
"^^{x.y) = v^(T^e"'^/(r), V2,±ie^^/(r)), where /(r) 
vanishes at r = and becomes 1 beyond a healing length 
^. This describes a ferromagnet (^ie"*"^, a/2, ±ze*'^)/2 
(at large distance) with a polar core (0, 1,0) (at r = 0). 
All spins line up in circles in the x?/-plane, with a mag- 
nitude |m| shrinks from 1 to zero as r ^ 0, (a "meron"). 
This is the state found by Yi and Pu.[5l. Fig. (la) and 
(lb) show the spin textures of ^ with upper and lower 
sign, which have opposite spin orientations. Note that re- 
versing spin orientation does not alter the phase winding 
(or circulation) of each spin component. 
(a) Hyperbolic planar texture: This is the same as (i) 
except with —cj). See Fig.(lc) and (Id). 
(Hi) Elliptic Mermin-Ho texture: This corresponds to 
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where /3(r) is an increasing function of r starting with 
/3(0) = 0. As r increases from to a distance where 
13 = 7r/2, C changes from (1,0,0) to (1,^726^^,-6^^^), 
which has the same circular spin texture as {i) at large 
distance. It differs from (i) in that it is everywhere ferro- 
magnetic. Vortex singularities can be eliminated without 
damaging the ferromagnetic order. See Fig.(le). 
(iv) Hyperbolic Mermin-Ho texture: This is the same as 
(in) with (j) —(j). 

(B) Energetic considerations: We shall discuss the 
zero field case at zero temperature. Although 5 7^ in 
current experiments, we first discuss this case because of 
its fundamental importance. In zero magnetic field, the 
energy is = T+V+ZY+Vd, where T = / |^V|^^p 



is the kinetic energy, V = / y(x)n(x) is harmonic trap. 



l/(x) = ^M{ujIx^ + coly^ + u;' 
describes density repulsion (cq > 0) and ferromagnetic in- 
teraction (c2 < 0)[3J, Vd = J 5'i(x)Ai(x-x')S'^ (x') 
is the dipolar energy, where Dij{x.) = Sij/x^ — 3xiXj/x^^ 
= 7^, 7 = QfI^b^ Of = —1/2 is the magnetic 
^-factor, and /i^ is the Bohr magneton. Noting that 

Ai(x) 



^(5ij(5(x), we have 



9d_ 
2 



/, 

y x,x' 



Q(x)Q(xO 

|x-x'| 



47r 

T 



L 



S^(x) 



(2) 



where Q(x) = V • S(x). The first term in eq.Q is pos- 
itive definite, the optimum spin configuration is the one 
that satisfies V • S = while keeping |m| = 1 to gain 
maximum ferromagnetic energy. The difficulty in finding 
the equilibrium textures is that these two conditions are 
not always compatible. 

For ^^Rb, C2/C0 = 0.005 pi, the healing length for spin 
is much longer than that for density, and exceeds the 
thickness of the "pancake" condensates. The spin degrees 
of freedom is effectively 2D, in the sense that 



-^^{x.y.z) = y^n{x,y,z)C^{x,y), 



(3) 



and m depends only on (x, y). Since n(x) is mirror sym- 
metric about the x?/-plane , the term QQ in eq.Q can 
be replaced by (5//(x)(5//(x') + Q_l(x)Qx(x'), where 
(5//(x) = mz{:sL±)dzn{x.), and (5±(x) = V± • Sx(x), 
V± = {Vx.Vy), S± = (Sx.Sy). The energy is mini- 
mized by rriz = Q± = 0, or simply -8 = with S in 
the xy-pldLiie. In regions where density is uniform, these 
conditions can be satisfied by the elliptic planar texture 
(i) mentioned above. Indeed, circular spin alignment 
are prevalent in all cases we studied. Their presence, 
however, often require some fractions of their hyperbolic 
counterparts to facilitate their close packing, even though 
the latter are not as energetically favorable. 

(C) Our calculation: For simplicity, we take the 
density to be a Gaussian along z and Thomas-Fermi 
in x?/-plane, n{x^y^z) = w{z)'^nFT{x^y)^ where w{z) = 
/2d ^1^^1/4^1/2^^ ^ -g width of the condensate 
along z, and nTF(r) = - ^(Mcj^x^ + Mcj^?/^)]/co, 
Co = Co/{V^d)^ r = (x,?/), and /i is the chemical 
potential determined by the total number of particles. 
We have ignored the effect of 62 on the density since 
C2/C0 = 0.005[l]. With this density, all the variational 
variables are contained in which can be parametrized 
by five real fields, {X^(r),z = 1..5}. To locate the en- 
ergy minima, we evolve the variables Xi{r) by a dissipa- 



, where F^ > 0. 



tive dynamics (iX^(r, t) /(it = ±. ^ ^x {y t) 
Since dE/dt = - J drri{SE/SXi{r,t))'^ < 0, this evolu- 
tion forces the energy to decrease, coming to a stop only 
when a local minimum {6E/SXi = 0) is reached. Start- 
ing this evolution with different initial conditions, one 
can locate the energy minimum while mapping out the 
energy surface in the space of Ci'^)^ o^" 
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FIG. 2: The middle section of the texture of an energy min- 
imum of a system with aspect ratio A = 10 in zero magnetic 
field. The state consists of a row of elliptic planar textures 
with alternating spin orientations. All spins lies in the xy- 
plane. See Fig. (la) and (lb). 

(D) Zero field case: We have performed imaginary 
time evolution for a system with 10^ particles for different 
aspect ratios A and different initial states: random config- 
urations, spiral spin textures, and vortex lattice textures. 
In all cases, the final state is either a uniform texture, or 
a state ^^^^ with a row of circular planar texture (type 
(i) in Section (A)) with alternating orientations. The 
texture for the A = 10 case is shown fig. 2. All the spins 
indeed lie in the x?/-plane as discussed in Section (B). 

Based on the discussions in (A), we find that our nu- 
merical result ^^^^ can be described accurately by 

y, z)-^ = V^(e'*(^)/(r), V2, e-^*(^)/(r)) (4) 

where n = w'^{z)nTF{^)i r = {x^y)^ ^(r) is the phase 
angle of the ID vortex lattice, 

^{x^y) = 7r/2 -\- diTgV{z)^ z = x-\-iy^ (5) 

and V{z) = Yl^{z - a - nb), n = 0, ±1, ±2, .., ±M is 
a polynomial describing an array of vortices of identical 
circulation centered at x = a, separated by distance b 
along the x axis. We find that when A is an odd (even) 
integer, a = (a = 6/2). M is the smallest integer such 
that {M-\-l)b > Rx- The function / describes the vortex 
core with size ^[7 . For all aspect ratios we examined, the 
optimal value of b (the size of the circular planar unit) 
and ^ (core size ) is found to be ~ 1.5 and ~ 0.25 in units 
of Ry (the shorter Thomas-Fermi in the xy-plane. 

The packing of vortices with identical vorticity in both 
ifji and if) -I may seem surprising, for it costs kinetic en- 
ergy. The reason is that this is the only packing that 
leads to circular spin textures, which is strong favored 
by dipolar energy. With dipolar energy scales as and 
kinetic energy scales as n, the condensate wavefunction 
will be determined by the former. Our calculations show 
that the vortex row state ^^^^ is the ground state for 
A < 2. For A > 2, the ground state is a uniform spin 
texture along x, while ^^^^ is a local minimum [8 . 

(E) A large magnetic field B: In a magnetic field 
B, the energy acquires linear and quadratic Zeeman shift, 

= -/7B ■ *;F^.*. and S2 = /</*;(B ■ ¥)%^,. 



For sufficiently large magnetic field, the spin texture of 
the system will rotate about B. In the frame rotating 
about B with Lamour frequency, the linear Zeeman shift 
is transformed away. Dipolar energy, however, acquires 
time dependent terms in this frame that oscillates with 
the Lamour frequency. Over times much longer than the 
Lamour cycle, the system sees an averaged dipolar energy 
% = /^,(r)A,,(r - v')Sj{v'), where 

_ (3B.B,-J.,) / _ 3(B.x)^ \ 
^mW- 2|x|3 y |xP )' 

The total time averaged energy is then E = T ^ 
Vd +^2- Note that E conserves B • / S. In ref.[l^ 2 , B is 
aligned with x up to a few degrees. The observed state 
with periodic spin texture is found to have j Sx ^ ^ and 
with most spins lie perpendicular to x. 

To search for the stationary states of E with J ~ 0, 
we performed the dissipative dynamics mentioned in Sec- 
tion (C) with a great variety of initial conditions. De- 
pending on the quadratic Zeeman energy (g), two types 
of stationary state emerge. For q/TiuOy > {q/%uOy)c ^ 0.1, 
the equilibrium state is is an uniaxial spin density wave 
with S(r) = zcos(i^B • r). For the parameters used in 
footnote [9], we find a wavelength 2t: / K ~ 12. bay ~ 
25/im, {ay = ^Jf^/ {MuOy)) and is independent of the an- 
gle between B and x. For q/fiuOy < {q/hujy)c, the system 
will settle in a state consisting of two large uniform spin 
domains along x and — x, with a pair of elliptic and hy- 
perbolic Mermin-Ho texture sandwiched in between. 

In additional to these two equilibrium states, we also 
find a class of spin textures which are essentially distorted 
vortex lattices which remain in the imaginary time evo- 
lution for a very long time, reflecting the flatness of the 
energy surface in the neighborhood of these states. The 
underlying structure of these distorted lattice are pairs 
of elliptic and hyperbolic MH vortices, (See Fig. 4), with 
essentially zero spin projection along x, J Sx/N ^ 10~^. 
These states remain long lived for q/Tiujy < 0.2 and be- 
gin to evolve more rapidly toward the stripe phase as q 
increases. We have calculated the spin-spin correlation 
function G(r) as defined in ref.[T for the vortex arrays 
in Fig. 4. The result is shown in Fig. 5 where bright 
color represents high value of G(r). The underlying lat- 
tice structure of the disordered spin texture Fig. 4 shows 
up a lattice of bright spots in Fig. 5. 

In Figure 6, we have plotted |6'(K)p, where S(K) 
is the Fourier transform of spin texture S(r) shown in 
Figure 4. The almost periodic structure in real space 
shows up as intense spots in iiT— space. Taking the 
brightest spot closest to the origin on the right, we find 

{ciy/^x-,ciy/\) ^ (0,0.08), or ^A^ + A^ ~ 25/im, com- 
parable to the period found in the stripe phase, and is the 
same order of the "dipolar length" [9 . We note that the 
corresponding plot in ref.pl shows a length scale ~ 10/im, 



FIG. 3: The middle section of the texture of the equihbrium 
state of E for a system with aspect ratio A = 10. 




FIG. 4: A long lived state of the time average energy ^ - an 
array of elliptic and hyperbolic Mermin-Ho vortex pair. 

differing from our result by a factor of 2. This difference 
may be due to our Gaussian approximation of the ac- 
tual density normal to the plane, which will contribute 
to systematic errors in the energies of all textures. How- 
ever, such systematic error will not change the fact that 
both the stripe phase and the spin lattice textures have 
similar length scales - a property that can be verified by 
experiments. 

Since the long lifetime in this dissipative evolution is 
due to the almost vanishing SE/S(^ similar long lifetime 
will also occur in real time evolution as it is driven by the 
same derivative. Our findings, which is consistent with 
those in ref.([Tl [2]), suggest that the observed periodic 
structure is a long lived dynamical state. They also pre- 
dict two equilibrium configurations: a row of alternating 
circular spin texture in zero field, and a single wavevector 
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FIG. 5: spin-spin correlation of the state in Fig. 4. Note that 
the bright spots form an almost square lattice. 
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FIG. 6: |S(K)|2 of the state in Fig.4; K = 27r(A-\A-^). 

spin stripe phase in large magnetic field. Our method is 
also applicable to other dipolar condensates . 
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